par(new=TRUE)
plot(WWNEnormal3$ASSET, WWNEnormal3$`distWWNE + distWWE`, type="l", col="red", lwd=1.5, xlim = c(0,100), ylim=c(0,0.005), lty=1, ylab="", xlab="")
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mar=c(4,4,2,2))
plot(JNEnormal1$ASSET, JNEnormal1$distJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00065), ylab="density", xlab="Wealth")
par(new=TRUE)
plot(JNEnormal2$ASSET, JNEnormal2$distJNE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00065), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
plot(JNEnormal3$ASSET, JNEnormal3$distJNE, type="l", col="red", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00065), lty=1, ylab="", xlab="")
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
plot(JNEnormal1$ASSET, JNEnormal1$distJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00025), ylab="density", xlab="Wealth")
par(mar=c(4,4,2,2))
plot(JNEnormal1$ASSET, JNEnormal1$distJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00025), ylab="density", xlab="Wealth")
par(new=TRUE)
plot(JNEnormal2$ASSET, JNEnormal2$distJNE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00025), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
plot(JNEnormal3$ASSET, JNEnormal3$distJNE, type="l", col="red", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00025), lty=1, ylab="", xlab="")
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
JNEnormal = read.table("distribution/distJ.out", header=TRUE)
JNEnormal1 = JNEnormal[JNEnormal$THETA == 1,]
JNEnormal1 = aggregate(distJNE ~ ASSET, JNEnormal1, sum)
JNEnormal2 = JNEnormal[JNEnormal$THETA == 2,]
JNEnormal2 = aggregate(distJNE ~ ASSET, JNEnormal2, sum)
JNEnormal3 = JNEnormal[JNEnormal$THETA == 0,]
JNEnormal3 = aggregate(distJNE ~ ASSET, JNEnormal3, sum)
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(WWNEnormal1$ASSET, WWNEnormal1$`distWWNE + distWWE`, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,0.005), ylab="density", xlab="Wealth")
par(new=TRUE)
plot(WWNEnormal2$ASSET, WWNEnormal2$`distWWNE + distWWE`, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,0.005), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
plot(WWNEnormal3$ASSET, WWNEnormal3$`distWWNE + distWWE`, type="l", col="red", lwd=1.5, xlim = c(0,100), ylim=c(0,0.005), lty=1, ylab="", xlab="")
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mar=c(4,4,2,2))
plot(JNEnormal1$ASSET, JNEnormal1$distJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00025), ylab="density", xlab="Wealth")
par(new=TRUE)
plot(JNEnormal2$ASSET, JNEnormal2$distJNE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00025), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
plot(JNEnormal3$ASSET, JNEnormal3$distJNE, type="l", col="red", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00025), lty=1, ylab="", xlab="")
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(WWNEnormal1$ASSET, WWNEnormal1$`distWWNE + distWWE`, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,0.005), ylab="density", xlab="Wealth")
par(new=TRUE)
plot(WWNEnormal2$ASSET, WWNEnormal2$`distWWNE + distWWE`, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,0.005), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
plot(WWNEnormal3$ASSET, WWNEnormal3$`distWWNE + distWWE`, type="l", col="red", lwd=1.5, xlim = c(0,100), ylim=c(0,0.005), lty=1, ylab="", xlab="")
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mar=c(4,4,2,2))
plot(JNEnormal1$ASSET, JNEnormal1$distJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00065), ylab="density", xlab="Wealth")
par(new=TRUE)
plot(JNEnormal2$ASSET, JNEnormal2$distJNE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00065), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
plot(JNEnormal3$ASSET, JNEnormal3$distJNE, type="l", col="red", lwd=1.5, xlim = c(0,100), ylim=c(0,0.00065), lty=1, ylab="", xlab="")
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mfrow=c(1,1))
par(mar=c(4,4,2,2))
plot(JNEnormal1$ASSET, JNEnormal1$distJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), ylab="density", xlab="Wealth")
par(new=TRUE)
plot(JNEnormal2$ASSET, JNEnormal2$distJNE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
plot(JNEnormal3$ASSET, JNEnormal3$distJNE, type="l", col="green", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), lty=1, ylab="", xlab="")
par(new=TRUE)
plot(JNEnormal4$ASSET, JNEnormal4$distJNE, type="l", col="orange", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), lty=1, ylab="", xlab="")
par(new=TRUE)
plot(JNEnormal5$ASSET, JNEnormal5$distJNE, type="l", col="pink", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), lty=1, ylab="", xlab="")
par(new=TRUE)
plot(JNEnormal6$ASSET, JNEnormal6$distJNE, type="l", col="red", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), lty=1, ylab="", xlab="")
par(new=TRUE)
plot(JNEnormal7$ASSET, JNEnormal7$distJNE, type="l", col="purple", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), lty=1, ylab="", xlab="")
abline(h=0, col="grey4", lty=2)
legend("topright", c("z=1", "z=2", "z=3", "z=4", "z=5", "z=6", "z=7"), col=c("black", "blue","green", "orange", "pink", "red", "purple"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
JNEnormal = read.table("distribution/distJ.out", header=TRUE)
JNEnormal1 = JNEnormal[JNEnormal$ZETA == 0,]
JNEnormal1 = aggregate(distJNE ~ ASSET, JNEnormal1, sum)
JNEnormal2 = JNEnormal[JNEnormal$ZETA == 1,]
JNEnormal2 = aggregate(distJNE ~ ASSET, JNEnormal2, sum)
JNEnormal3 = JNEnormal[JNEnormal$ZETA == 2,]
JNEnormal3 = aggregate(distJNE ~ ASSET, JNEnormal3, sum)
JNEnormal4 = JNEnormal[JNEnormal$ZETA == 3,]
JNEnormal4 = aggregate(distJNE ~ ASSET, JNEnormal4, sum)
JNEnormal5 = JNEnormal[JNEnormal$ZETA == 4,]
JNEnormal5 = aggregate(distJNE ~ ASSET, JNEnormal5, sum)
JNEnormal6 = JNEnormal[JNEnormal$ZETA == 5,]
JNEnormal6 = aggregate(distJNE ~ ASSET, JNEnormal6, sum)
JNEnormal7 = JNEnormal[JNEnormal$ZETA == 6,]
JNEnormal7 = aggregate(distJNE ~ ASSET, JNEnormal7, sum)
par(mfrow=c(1,1))
par(mar=c(4,4,2,2))
plot(JNEnormal1$ASSET, JNEnormal1$distJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), ylab="density", xlab="Wealth")
par(new=TRUE)
plot(JNEnormal2$ASSET, JNEnormal2$distJNE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
plot(JNEnormal3$ASSET, JNEnormal3$distJNE, type="l", col="green", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), lty=1, ylab="", xlab="")
par(new=TRUE)
plot(JNEnormal4$ASSET, JNEnormal4$distJNE, type="l", col="orange", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), lty=1, ylab="", xlab="")
par(new=TRUE)
plot(JNEnormal5$ASSET, JNEnormal5$distJNE, type="l", col="pink", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), lty=1, ylab="", xlab="")
par(new=TRUE)
plot(JNEnormal6$ASSET, JNEnormal6$distJNE, type="l", col="red", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), lty=1, ylab="", xlab="")
par(new=TRUE)
plot(JNEnormal7$ASSET, JNEnormal7$distJNE, type="l", col="purple", lwd=1.5, xlim = c(0,100), ylim=c(0,0.0002), lty=1, ylab="", xlab="")
abline(h=0, col="grey4", lty=2)
legend("topright", c("z=1", "z=2", "z=3", "z=4", "z=5", "z=6", "z=7"), col=c("black", "blue","green", "orange", "pink", "red", "purple"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
CollateralJNE = read.table("VF/JNEvalue.out", header=TRUE)
CollateralJNE = read.table("VF/JNEvalues.out", header=TRUE)
CollateralJNE
CollateralJE = read.table("VF/JEvalues.out", header=TRUE)
CollateralJNE = CollateralJNE[CollateralJNE$THETA == 1 && CollateralJNE$Z == 4,]
par(mfrow=c(1,1))
par(mar=c(4,4,2,2))
plot(CollateralJNE$ASSET, CollateralJNE$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,1.5), ylab="Search effort (sw)", xlab="Wealth")
CollateralJNE
plot(CollateralJNE$ASSET, CollateralJNE$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,1.5), ylab="Search effort (sw)", xlab="Wealth")
CollateralJNE$collateralJNE
CollateralJNE$THETA == 1
CollateralJNE = read.table("VF/JNEvalues.out", header=TRUE)
CollateralJE = read.table("VF/JEvalues.out", header=TRUE)
CollateralJNE = CollateralJNE[CollateralJNE$THETA == 1 && CollateralJNE$Z == 4,]
CollateralJNE
CollateralJNE = read.table("VF/JNEvalues.out", header=TRUE)
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 1 && CollateralJNE$Z == 4,]
CollateralJNE$THETA == 1
CollateralJNE$Z == 4
CollateralJNE$THETA == 1 && CollateralJNE$Z == 4
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 1 & CollateralJNE$Z == 4,]
CollateralJNE1
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,1.5), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,1000), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylab="Search effort (sw)", xlab="Wealth")
CollateralJE1 = CollateralJE[CollateralJE$THETA == 1 & CollateralJE$Z == 4,]
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,1.5), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylim=c(0,100), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylim=c(0,150), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,150), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylim=c(0,150), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,150), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
plot(UiNEnormal3$ASSET, UiNEnormal3$searchUiNE_W, type="l", col="red", lwd=1.5, xlim = c(0,100), ylim=c(0,1.5), lty=1, ylab="", xlab="")
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylim=c(0,150), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,150), ylab="", xlab="", axes=FALSE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,300), ylim=c(0,150), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylim=c(0,150), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,300), ylim=c(0,150), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,150), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,150), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylim=c(0,150), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,150), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 1 & CollateralJNE$Z == 6,]
CollateralJE1 = CollateralJE[CollateralJE$THETA == 1 & CollateralJE$Z == 6,]
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,100), ylim=c(0,150), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,150), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylim=c(0,150), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylim=c(0,300), ylab="Search effort (sw)", xlab="Wealth")
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 1 & CollateralJNE$Z == 5,]
CollateralJE1 = CollateralJE[CollateralJE$THETA == 1 & CollateralJE$Z == 5,]
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylim=c(0,300), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,150), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylim=c(0,500), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,400), ylim=c(0,500), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,400), ylim=c(0,500), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,400), ylim=c(0,150), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,400), ylim=c(0,400), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,400), ylim=c(0,400), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,250), ylim=c(0,400), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,250), ylim=c(0,400), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylim=c(0,400), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,300), ylim=c(0,400), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,350), ylim=c(0,400), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,350), ylim=c(0,400), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,300), ylim=c(0,150), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,100), ylim=c(0,150), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,5), ylim=c(0,600), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,20), ylim=c(0,600), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,50), ylim=c(0,60), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,20), ylim=c(0,60), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,20), ylim=c(0,20), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,10), ylim=c(0,10), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,10), ylim=c(0,10), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,10), ylim=c(0,10), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,3), ylim=c(0,3), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,3), ylim=c(0,3), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,3), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,2), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,3), ylim=c(0,3), ylab="", xlab="", axes=FALSE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,3), ylim=c(1.0,1.5), ylab="", xlab="", axes=FALSE)
7par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,3), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,3), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
legend("topright", c("int.rate","loan"), col=c("blue", "black"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("int.rate","loan"), col=c("blue", "black"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(0.9,1.7), ylab="", xlab="", axes=FALSE)
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(0.9,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("int. rate","loan"), col=c("blue", "black"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 1 & CollateralJNE$Z == 3,]
CollateralJE1 = CollateralJE[CollateralJE$THETA == 1 & CollateralJE$Z == 3,]
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,350), ylim=c(0,400), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,350), ylim=c(0,400), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("int. rate","loan"), col=c("blue", "black"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,5.7), ylab="", xlab="", axes=FALSE)
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 1 & CollateralJNE$Z == 5,]
CollateralJE1 = CollateralJE[CollateralJE$THETA == 1 & CollateralJE$Z == 5,]
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,350), ylim=c(0,400), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,350), ylim=c(0,400), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("int. rate","loan"), col=c("blue", "black"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
CollateralJNE1$rateJNE
abline(h=0, col="grey4", lty=2)
abline(v=1.2, col="grey4", lty=2)
axis(side = 4)
mtext(side = 4, line = 3, 'Number genes selected')
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 0 & CollateralJNE$Z == 5,]
CollateralJE1 = CollateralJE[CollateralJE$THETA == 0 & CollateralJE$Z == 5,]
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,350), ylim=c(0,400), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,350), ylim=c(0,400), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(v=1.2, col="grey4", lty=2)
axis(side = 4)
mtext(side = 4, line = 3, 'Interest rate')
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="loan", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(v=1.2, col="grey4", lty=2)
axis(side = 4)
mtext(side = 4, line = 3, 'interest rate')
legend("topright", c("int. rate","loan"), col=c("blue", "black"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 2 & CollateralJNE$Z == 5,]
CollateralJE1 = CollateralJE[CollateralJE$THETA == 2 & CollateralJE$Z == 5,]
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,350), ylim=c(0,400), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,350), ylim=c(0,400), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="loan", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(v=1.2, col="grey4", lty=2)
axis(side = 4)
mtext((1,1), line = 3, 'interest rate')
text((1,1), line = 3, 'interest rate')
text((1), line = 3, 'interest rate')
text((1),  'interest rate')
text((1),1,  'interest rate')
text(1,1,  'interest rate')
text(x=1,y=1,'interest rate')
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="loan", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(v=1.2, col="grey4", lty=2)
text(x=1,y=1,'interest rate')
axis(side = 4)
mtext(side = 4, line = 3, 'interest rate')
legend("topright", c("int. rate","loan"), col=c("blue", "black"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
text(x=1,y=2,'interest rate')
text(x=1,y=1.5,'interest rate')
text(x=0.5,y=1.5,'interest rate')
text(x=0.5,y=1.5,'exemption level')
text(x=0.65,y=1.5,'exemption level')
text(x=0.65,y=1.7,'exemption level')
text(x=0.65,y=1.6,'exemption level')
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="loan", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(v=1.2, col="grey4", lty=2)
text(x=0.65,y=1.6,'exemption level')
axis(side = 4)
mtext(side = 4, line = 3, 'interest rate')
legend("topright", c("int. rate","loan"), col=c("blue", "black"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
CollateralJNE = read.table("VF/JNEvalues.out", header=TRUE)
CollateralJE = read.table("VF/JEvalues.out", header=TRUE)
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 2 & CollateralJNE$Z == 5,]
CollateralJE1 = CollateralJE[CollateralJE$THETA == 2 & CollateralJE$Z == 5,]
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,350), ylim=c(0,400), ylab="Search effort (sw)", xlab="Wealth")
## POLICY FUNCTION :: COLLATERAL VALUE ##
CollateralJNE = read.table("VF/JNEvalues.out", header=TRUE)
CollateralJE = read.table("VF/JEvalues.out", header=TRUE)
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 2 & CollateralJNE$Z == 5,]
CollateralJE1 = CollateralJE[CollateralJE$THETA == 2 & CollateralJE$Z == 5,]
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,350), ylim=c(0,2), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,350), ylim=c(0,2), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,2), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,3), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,5), ylab="Search effort (sw)", xlab="Wealth")
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,4), ylab="Search effort (sw)", xlab="Wealth")
CollateralJNE = read.table("VF/JNEvalues.out", header=TRUE)
CollateralJE = read.table("VF/JEvalues.out", header=TRUE)
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 2 & CollateralJNE$Z == 5,]
CollateralJE1 = CollateralJE[CollateralJE$THETA == 2 & CollateralJE$Z == 5,]
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,4), ylab="Search effort (sw)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(0,4), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("Low-skilled","Middle-skilled", "High-skilled"), col=c("red", "black","blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
legend("topright", c("Access to credit","Excluded"), col=c("black", "blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,4), ylab="capital invested (k)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(0,4), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("Access to credit","Excluded"), col=c("black", "blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="loan", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(v=1.2, col="grey4", lty=2)
text(x=0.65,y=1.6,'exemption level')
axis(side = 4)
mtext(side = 4, line = 3, 'interest rate')
legend("topright", c("int. rate","loan"), col=c("blue", "black"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
legend("topright", c("int. rate","loan"), col=c("black", "black"), lwd=1.5, cex=0.8, lty=c(1,5,1), y.intersp=1.2)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", lty = 5, col="black", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="loan", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", lty = 5, col="black", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(v=1.2, col="grey4", lty=2)
text(x=0.65,y=1.6,'exemption level')
axis(side = 4)
mtext(side = 4, line = 3, 'interest rate')
legend("topright", c("int. rate","loan"), col=c("black", "black"), lwd=1.5, cex=0.8, lty=c(1,5,1), y.intersp=1.2)
## POLICY FUNCTION :: COLLATERAL VALUE ##
CollateralJNE = read.table("VF/JNEvalues.out", header=TRUE)
CollateralJE = read.table("VF/JEvalues.out", header=TRUE)
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 2 & CollateralJNE$Z == 5,]
CollateralJE1 = CollateralJE[CollateralJE$THETA == 2 & CollateralJE$Z == 5,]
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,4), ylab="capital invested (k)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(0,4), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("Access to credit","Excluded"), col=c("black", "blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="loan", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", lty = 5, col="black", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(v=1.2, col="grey4", lty=2)
text(x=0.65,y=1.6,'exemption level')
axis(side = 4)
mtext(side = 4, line = 3, 'interest rate')
legend("topright", c("int. rate","loan"), col=c("black", "black"), lwd=1.5, cex=0.8, lty=c(1,5,1), y.intersp=1.2)
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="loan", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", lty = 5, col="black", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(v=1.2, col="grey4", lty=2)
text(x=0.75,y=1.6,'exemption level')
axis(side = 4)
mtext(side = 4, line = 3, 'interest rate')
legend("topright", c("int. rate","loan"), col=c("black", "black"), lwd=1.5, cex=0.8, lty=c(1,5,1), y.intersp=1.2)
## POLICY FUNCTION :: COLLATERAL VALUE ##
CollateralJNE = read.table("VF/JNEvalues.out", header=TRUE)
CollateralJE = read.table("VF/JEvalues.out", header=TRUE)
CollateralJNE1 = CollateralJNE[CollateralJNE$THETA == 2 & CollateralJNE$Z == 5,]
CollateralJE1 = CollateralJE[CollateralJE$THETA == 2 & CollateralJE$Z == 5,]
par(mfrow=c(1,2))
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$collateralJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,4), ylab="capital invested (k)", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJE1$collateralJE, type="l", col="blue", lwd=1.5, xlim = c(0,2), ylim=c(0,4), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(h=0, col="grey4", lty=2)
legend("topright", c("Access to credit","Excluded"), col=c("black", "blue"), lwd=1.5, cex=0.8, lty=c(1,1,1), y.intersp=1.2)
par(mar=c(4,4,2,2))
plot(CollateralJNE1$ASSET, CollateralJNE1$SloanJNE, type="l", col="black", lwd=1.5, xlim = c(0,2), ylim=c(0,1.5), ylab="loan", xlab="Wealth")
par(new=TRUE)
plot(CollateralJE1$ASSET, CollateralJNE1$rateJNE, type="l", lty = 5, col="black", lwd=1.5, xlim = c(0,2), ylim=c(1.0,1.7), ylab="", xlab="", axes=FALSE)
par(new=TRUE)
abline(v=1.2, col="grey4", lty=2)
text(x=0.75,y=1.6,'exemption level')
axis(side = 4)
mtext(side = 4, line = 3, 'interest rate')
legend("topright", c("int. rate","loan"), col=c("black", "black"), lwd=1.5, cex=0.8, lty=c(1,5,1), y.intersp=1.2)
